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0^ ■ Abstract. We show that the diagonal matrix elements (Or p ), where O = {1, (3, iocn/3} are the 



o 
o 

(N 



standard Dirac matrix operators and the angular brackets denote the quantum-mechanical average 
for the relativistic Coulomb problem, may be considered as difference analogs of the radial wave 
functions. Such structure provides an independent way of obtaining closed forms of these matrix 
^ ' elements by elementary methods of the theory of difference equations without explicit evaluation of 

the integrals. Three-term recurrence relations for each of these expectation values are derived as a 
by-product. Transformation formulas for the corresponding generalized hypergeometric series are 
discussed. 
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^ ' 1. Introduction 

a 

Recent experimental and theoretical progress has renewed interest in quantum electrodynamics 
of atomic hydrogenlike systems (see, for example, [8], [9], p2], [Ej, ^5\, [22], and [21] and references 
therein). In the last decade, the two-time Green's function method of deriving formal expressions for 
the energy shift of a bound-state level of high-Z few-electron systems was developed [22] and numer- 
ical calculations of QED effects in heavy ions were performed with excellent agreement to current 
experimental data [8], [9], [21]. These advances motivate detailed study of the expectation values of 
the Dirac matrix operators between the bound-state relativistic Coulomb wave functions. Special 
cases appear in calculations of the magnetic dipole hyperfine splitting, the electric quadrupole hy- 
perfine splitting, the anomalous Zeeman effect, and the relativistic recoil corrections in hydrogenlike 
ions (see, for example, p], |21j . [23] . [25] and references therein). These expectation values can be 
used in calculations with hydrogenlike wave functions when a high precision is required. 

X. 

In the previous paper [25] . we have evaluated the relativistic Coulomb integrals of the radial 



O 



functions, 



POO 

A p = / r p+2 (F 2 (r) + G 2 (r)) dr, (1.1) 
Jo 

POO 

B p = r p+2 [F 2 (r) - G 2 (r)) dr, (1.2) 



DC 



C p = / r p+2 F{r)G{r) dr, (1.3) 



o 



for all admissible powers p, in terms of three special generalized hypergeometric 3F2 series related to 
the Chebyshev polynomials of a discrete variable [16J (we concentrate on the radial integrals since, 
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for problems involving spherical symmetry, one can reduce all expectation values to radial integrals 
by use of the properties of angular momentum). These integrals are linearly dependent: 

(2k + e(p+ 1)) A p - (2ek + p + 1) B p = 4pC p (1.4) 

(see, for example, [I], [19], [2D] , and [25] for more details). Thus, eliminating, say C p , one can deal 
with A p and B p only. The corresponding representations in terms of only two linearly independent 
generalized hypergeometric series are given in this paper (see (j3.1j) - (|3.3p and (j3.4j) - (l3.6p ). 

The integrals (ll.ll) - fll.3p satisfy numerous recurrence relations in p, which provide an effective 
way of their evaluation for small p (see [T], [19], [20], [25] and references therein). The two-term 
recurrence relations were derived by Shabaev [19], [20] on the basis of a hypervirial theorem and 
by a different method using relativistic versions of the Kramers-Pasternack three-term recurrence 
relations in [26 1. In our notations, 



Au 2 e + 2K(p + 2) + e(p+l)(2K6+p + 2) 
Ap+1 ~ ~ {P + 1) 4(1 -e 2 ) (p + 2)#z Ap (L5) 

4p 2 (p + 2) + (p + 1) (2k£ + p + l) [2k£ + p + 2) 
+ A(l-e 2 )(p + 2)(3iu p ' 



and 



o , , 1 , 4i/ 2 + 2Kg(2p + 3)+e 2 (p+l)(p + 2) 

4/A (p + 2) + (p + 1) (2/ce + p + 1) (2k + 6 (p + 2)) 
+ 4(l-e 2 )(p + 2)/3p p 

Vg (p + 1) + p (2/t £ + p) (2k + e(p + 1)) 
fi (Aw 2 — p z )p 
4/i 2 (p + 1) + p (2ke + p) (2ke + p + 1) 
jj, [Av z — p z ) p 

n Av 2 + 2KS (2p + l)+ £ 2 p(p+l) 
fi (4u z — p z ) 



} Au 2 e + 2k (p + 1) + ep (2ke + p + 1) 
p. (4i/ 2 — p 2 



~TJ~~~9 3\ 



respectively. Here, 



« = ± (j + 1/2) , !/ = ^K 2 -fI 2 , 



fi = aZ = Ze 2 /hc, a = yl — e 2 , (1.9) 
e = E/mc 2 , (3 = mc/h 
with the total angular momentum j = 1/2,3/2,5/2, ... (see |25j and [27] for more details). 

These recurrence relations are complemented by the symmetries of the integrals A p , B p , and C p 
under reflections p — ► — p — 1 and p — > — p — 3 found in [25] (see also [2]). For example, 

A_ p _ 3 = (2a(3) w 1 {2V ' P ~ 2 \ (1.10) 

p v p; r(2i/ + p + 3) v 7 
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x f_ (p + 1) ^ + 2 £ K(2p + 3)-(p + 2f 
^ ^ ; p+2 ? 



2/t (2£K - 1) P+ B„ . . 

p + 2 p 1 ' 



o p+3 r(2z/-p-2) 



E_ p _ 3 = (2^)^ ^ ^— =( (1.11) 

p v M; T(2^ + p + 3) v ; 

x (e (p + 1) (2p + 3) A p + (4z/ 2 - 2sk (2p + 3) - (p + l) 2 ) 5 P ) , 
for independent convergent integrals A p and B p . 

In this paper, we would like to draw reader's attention to an interesting analogy between the 
explicit solutions of the first order system of difference equations fll.5ft - fll.6p and the standard 
method of dealing with the system of differential equations for the radial relativistic Coulomb wave 
functions F and G (see, for example, [1], [5], [10], [17], [18], and [27J regarding solution of the Dirac 
equation in Coulomb field). En route, we derive the three-term recurrence relations for each of the 
single integrals fll.ip - (jl.3ft that seem to be new and convenient for their evaluation. Our observation 
provides an independent method of obtaining closed forms of these matrix elements, but this time, 
from the theory of difference equations and without explicit evaluation of the integrals. Some 
transformation formulas for the corresponding generalized hypergeometric series are derived. 



2. Three-term Recurrence Relations 

Several relativistic Kramers-Pasternack three-term vector recurrence relations for the integrals 
-Ap^ Bp , Cp have been obtained in [26]. A more general setting is as follows. Let us rewrite f ll.5p - fll.6l) 
and fll.7p - fjl.8ft in the matrix form 

A P \ JAY t V, \ = t A \ (2a) 



and denote 



with 



Bp J \ B p _i J \ Bp^i J \ Bp 



Cp Up I ^p \ Cp a p 



4v 2 e + 2k (p + 1) + ep (2ks + p + 1) 
ap ~ P 4(1 -e 2 ){p+l)^ ' {2 - 6) 

V (p + 1) + p (2k€ + p) {2KB + P+1) 

4 (1 - e 2 ) (p + 1) ' K ] 

_ 4z/ 2 + 2ke (2p + 1) + e 2 p (p + 1) 
° P ~ ~ P 4(l- e 2)(p+i)^ ' W 

_ 4/i 2 e (p+l)+p {2ks + p) {2k + e(p + 1)) 
p ~ 4(l-£ 2 )(p+l)^ 



4 

and 
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A p = det S p 



(4i/ 2 -p 2 )p 



(2.7) 

(2a(3) 2 (p + l) 

Eliminating A p and B p , respectively, from the system (12. ip . we arrive at the following three-term 
recurrence equations for the independent integrals 



.4 



p+i 



Bp+i 



a p+l ~\ ^ a p j sip ^ ^p ^ip-l) 

dp+i H — ^— cip ) Bp ^—A p Bp-i, 



(2.8) 
(2.9) 



which seem are missing in the available literature. 



In general, one can easily verify that the following vector three-term recurrence relation holds: 



Ap+i 
B, 



Dr. 



for two matrices M p and N p provided that 

S p+1 = M P + NpS p \ 



Bp-i 



(2.10) 



(2.11) 

Our equations (I2.8p - fl2.9p provide a diagonal matrix solution. According to (12. ip . (12.101) and (12. lip , 
a simple identity 

' A r 



B,, 



(Sp+i — NpSp x ) 



(2.12) 



,+AU P ~ l 

} p+i J \ Bp J \ Bp_i 

holds for any matrix N p . The known three-term recurrence relations for the relativistic Coulomb 
integrals can be obtained by choosing different forms of the matrix N p . The case N p = goes back 
to the two-term recurrence relation (12.11) and two more explicit solutions have been found in 
Here we analyze another possibility and take 



X P 



and N„ 



A r 



Hp J \ Vp o 

for suitable parameters X p and A new convenient relations are as follows 



A 
B 



p+i 



p+i 







-A r 



' J p+i/b p 







, a p 

"p+1 ~r Cp+i — . 

Cp / 







A p 
B„ 



Gp+l/cp 



B 



(2.13) 



p-i 



A 
B 



p+i 



p+i 



o-p+i + b p+ i- 



ip+i + c p+i~r 



p_ 

dp ) 



A p 
B,. 
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A. 
B 



p+i 



p+i 



Jp+i + a p+i 



Op \ 



c p+i + dp+i 



dr. 



Ap 



-A r 



dp+i/bp 



a p+i/ c p 




A p -i 



B, 



p-i 



Explicit diagonal form, when the equations are separated, is given by 

HP(p) 



A 



p+i 



a? (5 (An 2 (p + l)+p (2ks + p) (2ks + p + 1)) (p + 2) 



A r . 



(Av 2 - p 2 ) (V (p + 2) + (p + 1) (2ks + p + 1) (2/ce + p + 2)) p 
(2a/?) 2 (p + 1) + p (2k £ + p) (2k£ + p + 1)) (p + 2) 



.4 



p— i) 



B, 



•p+i 



efxQ (p) 



Br, 



where 



a 2 (3 (Av 2 + 2ks (2p + 1) + £ 2 p (p + 1)) (p + 2) 
(Av 2 - p 2 ) (4z/ 2 + 2ke (2p + 3) + e 2 (p + 1) (p + 2)) (p + 1) 
(2a/?) 2 (4i/ 2 + 2k£ (2p + 1) + e 2 p (p + 1)) (p + 2) 

P (p) = 2ep A + (8ke 2 + he) p 3 

+ (8k 2 £ 3 + 8k 2 s + 20ks 2 - 4k - 8z/ 2 £ + 3e) p 2 
+ (12kV + 20n 2 e + 16ne 2 - 10k - 20v 2 e) p 
+An 2 E i + 8k 2 e + Ake 2 - 4k - 12z/ 2 e, 

Q (p) = 2e 2 p 3 + (7e 2 + 8ke - 2) p 2 
+ (8v 2 + 7e 2 + IQke - A) p 
+ 12v 2 + 2e 2 + Que - 2. 



B 



p—ii 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



In comparison with other papers (see pQ, [2], [TO], [20], [23], [20] , and references therein), our 
consideration provides an alternative way of the recursive evaluation of the special values A p and 
B p , when one deals separately with one of these integrals only. The corresponding initial data can 
be found in [23J. It is important emphasizing, for the purpose of this paper, that this argument 
resembles the reduction of the first order system of differential equations for relativistic radial 
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Coulomb wave functions F and G to the second order differential equations (see, for example, [17] 
and [27J). 

If one wants to solve equations fl2.17l) - fl2.18l) analytically for all admissible powers, then the 
major obstacle is that they are not difference equations of hypergeometric type on a quadratic 
lattice, solutions of which are available in the literature [16] . The following consideration helps. A 
linear transformation 

X ' ^ T p ( A 2 I , (2.21) 



where 



T P = ( Qp ^ P 1 , det T p = a p 5 p - (3 plp ? 0, (2.22) 



1 P 

results in a new system of the first order difference equations 



with 



where the corresponding similar matrix is given by 

S p = T p S p T p _\ =(l P M (2.24) 

\ Up j 

det T p _! a p = a p 5 p „ia p - a p 7 p _ 1 fe p + f3 p 5 p -xc p - fifi v _ x d p , (2.25) 

det T p _i b p = -a p {5 p _ x a p + a v a v -\b p - P> v fi v _ x c p + I3 p a p -id p , (2.26) 

det T p _i Bp = 7 p (5 p _ia p - 7 p 7 P -i & P + VWCp - 5 p ^ p ^ 1 d p , (2.27) 

det T p _! ^ = -jpPp^cip + -i p a p „ib p - S^^Cp + 5pa p _i<i p , (2.28) 

i j rji 

A p = det S p = det 5 P 6 p ^ 0. (2.29) 

det 

The new separated three-term recurrence equations take the similar forms 

X p+ i = ( a p+ i + ~ + rfp ) A p ^— Ap Xp_i, (2.30) 



and 




Yp+i — I d p+ i H — — — a p I Yp ^— Ap yp_i- (2.31) 



As in the case of the radial wave functions [T7] and [27], there are several possibilities to choose 
the matrix T p in order to simplify the original equations (I2.17I) - (I2.18I) . Examples of such transfor- 
mations, when the resulting equations are of a hypergeometric type and coincide with difference 
equations for special dual Hahn polynomials [TT], [2], [TB] (see also appendix A), are given in the 
next section. 
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3. Transformations of Relativistic Coulomb Integrals 

The integrals A p , B p , and C p can be evaluated in terms of two linearly independent 3F2 func- 
tions, which are related to the special dual Hahn polynomials that can be thought of as difference 
analogs of the Laguerre polynomials in explicit formulas for the radial wave functions (see [17] and 
[27] for a detailed tutorial on solution of the relativistic Coulomb problem). This fact has been 
partially explored in [22] and we elaborate on this connection here. Two different representations 
of the expectation values are available in a complete analogy with the well-known structure of the 
relativistic wave functions. 

Analogs of the traditional forms are as follows 

2 (p + 1) o/i (2aPf r ^ + ^ n A p (3.1) 

r (2z/ + p + 1) 

= (/i + an) (a (2ek + p + 1) — 2e/i) 
1 — n, —p, p + 1 
2z/ + l, 1 
+ (pi — a/c) (a (2sk + p + 1) + 2efx) 
-n, -p, p + 1 
2^ + 1, 1 



x 3^2 



x 3^2 



2 (p + 1) a/i (2a/3) p Wo \ , 3 5 P (3.2) 



x 



3^2 



r(2z/ + i) 
r(2z/ + P + i) 

(/i + a/t) (a (2k + e (p + 1)) — 2^) 
1 — n, — p, p + 1 
2z/ + l, 1 
+ (w - aft) (a (2k + £ (p + 1)) + 2/i) 
—n, —p, p + l 
2u + l, 1 

Ap{2a(3f ' ; } . C p (3.3) 





The averages of r p for the relativistic hydrogen atom were evaluated in the late 1930s by Davis [B] 
as a sum of certain three 3F2 functions. But it has been realized only recently that these series 
are, in fact, linearly dependent and related to the Chebyshev polynomials of a discrete variable 
[25] . Our formulas here present the final result (we use the standard definition of the generalized 
hypergeometric series throughout the paper [3], [7J). 
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Analogs of the Nikiforov and Uvarov form of the relativistic radial functions [17], [27] are given 

by 

4 (p + 1) e\w (2a(3f A p (3.4) 
= a (en + v) (2 (ek — i/)+p + l) 

T(2u + p + 3) ( 1 - n, p + 2, -p - 1 
X T(2z/ + 2) 3 2 ^ 2z/ + 2, 1 
—a (e« — u) (2 (e/c + z/) + p + 1) 

r (2v + p + 1) ( -n, p + 2, -p-1 

x r (2 ,) ^ 2,, 1 

4/iz/(2a/3) p 5 P (3.5) 

, . T(2z/ + p + 3) / l-Ti, p + 2, 
= a £K + z/ ; — 3 F 2 

V J T(2u + 2) 3 2 y 2z/ + 2, 1 

. T(2z/ + p+l) /-n, p + 2, -p-1 

— a [EK — V) — r 

V ; T (2z/) 3 2 y 2u, 1 

8 (p + 1) ^ 2 z/ (2a/5) p C p (3.6) 
= a (ek + v) (2k (ek — u) + (p + 1) (k — ez/)) 
r(2i/ + p + 3) / 1-n, p + 2, -p-1 
T(2z/ + 2) 3 2 ^ 2z/ + 2, 1 
—a (ek — v) (2k (ek + u) + (p + 1) (k + ei/)) 
T(2z/ + p + l) / -n, p + 2, -p-1 

X ; : iro 

r (2//) 3 2 { 2i/, 1 

These representations simplify Eqs. (3.7)-(3.9) of [25] with the help of the linear relation (jl.4p (the 
calculation details are left to the reader). 

It is important noting in this paper, that formulas (I3.ip - (l3.3p and (I3.4l) - (l3.6p provide explicit 
examples (of inverses) of the linear transformations (12.211) that reduce the original three-term recur- 
rence relations (I2.8p - (l2.9p to the difference equations of the corresponding dual Hahn polynomials 
in a complete analogy with the case of the relativistic radial wave functions (see, for example, [TTJ 
and [27J). One may choose any two of three linearly dependent integrals A p , B p , and C p and take 
the corresponding renormalized dual Hahn polynomials as X p and Y p . 

For example, by choosing A p and B p as the independent integrals and introducing 

X, = 3Fi 1 -"\7' P , + i y (37) 



Y, 



p — 3-^2 



2Z/+1, 1 

-n, -p, p + 1 
2z/+l, 1 
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from (I3.ip - fl3.2l) we arrive at the following transformation matrix 

(2a(3) p T (2v + 1) 



T 



2a 2 ro + p+i) 

( a (2k + e(p+l)) + 2/i a (2sk + p + 1) + 2e/x \ 

p + a« H + an 

a(2/c + e(p+l))-2/x a (2eK + p + 1) - 2e/x 



with 



Then 



V 



det T„ 



fi — an 



(2a(3) p 



r (2i/ + 1) 



/j, — aK 



li(p + l) 



J 



r(2z/ + p+l)y a(/i 2 -a 2 K 2 )' 



5*1, 



-i 



x 



T^T-_\= (a 2 p(2z/+p)) 
—a 2 p 2 + 2ae/jp — 2 (/j 2 
-2 (p 2 - a 2 K 2 ) 

with the help of the matrix identity (IB.ip and 

A p = det S p 

The new system (12.231) takes much simplier form 



2 (/x 2 - a 2 ^ 



a 2 p 2 + 2a£pp + 2 (fx 2 — <x 2 k 2 ) 



2u — p 
2v + p' 



a 2 p 2 — 2ae\xp + 2 (p 2 — a 2 K 2 ) 

a 2 p (2v + p) 
2 (/i 2 - a 2 /t 2 ) 



p-i 



and 



a 2 p(2z/ + p) 
2(p 2 



Yp-i 



a 2 K 2 \ 



a 2 p (2v + p) v 1 
a 2 p 2 + 2aepp + 2 (/i 2 — a 2 /t 2 ) 



a 2 p (2v + p) 



with the initial data 

Y 
1 

: 2d 2 



A 
Bo 



/ a(2K + e) + 2fi a (2ek + 1) + 2ep \ 

p + an p + an 

a(2K + e)-2/j, a (2sk + 1) - 2ep 

V p — <xk p — an / 



1 

e 



(3.8) 



(3.9) 
(3.10) 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



p — an p — an 

After this transformation, the three-term recurrence relations f l2.30p -f l2.3ip coincide with the dif- 
ference equations for the corresponding special dual Hahn polynomials ( I A. 71) (one should use the 
spectral identity ep = a(u + n), further computational details are left to the reader). Our con- 
sideration shows how the relativistic Coulomb expectation values A p and B p can be independently 
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obtained in their closed forms (I3.ip - fl3.2l) . when solving the original system fll.5l) - fll.6l) by the 
methods of the theory of difference equations developed in the previous section and without explicit 
evaluation of the integrals. A striking similarity with the structure of the radial wave functions 
provides a guidance in this approach. A similar analysis of the case f l3.4j) -f l3T5{) is left to the reader. 

On the second hand, our equations ( 13.1l) -f j3~3l and ( 13.4p -( l3T6l) imply the following linear relations: 

^ ( ±-n, -p, p+l\ 
^ 2v + l, 1 J (3 - 15) 
_ {2v + n){2v + p + 1) (2u + p + 2) (2w + p + 1) 
4i/ {2v + 1) (y + n) (p + 1) 
1 — n, p + 2, —p — 1 
2z/ + 2, 1 



and 




-n, p + 2, -p - 1 
2(i/ + n)(p+l) 3 2 { 2v, 1 



3 F 2 [ + 1 , (3.16) 

" 1 2u + 1, 1 1 

n (4v + 2n - p - 1) (2v + p + 1) {2v + p + 2) 



4z/ (2z/ + 1) (u + n) (p + 1) 
1 — n, p + 2, — p — 1 
2z/ + 2, 1 



[2i/ + n) (2n-p- 1) ( -n, p + 2, -p - 1 

3-^2 



2(i/ + n)(p+l) ^ 2i/, 1 

between two pairs of the generalized hypergeometric series under consideration. As required, only 
one dimensionless parameter is involved in the transformations. Details of these elementary but 
rather tedious calculations are left to the reader. 



In addition, from (3.7) of [25] and (13.41) of this paper one gets 



p (p + 1) / l-n, p + 1, -p 
2u + n 3 2 I 2i/ + l, 2 



(3.17) 



(p-2u) {2v + p+l) f 1-n, p + 1, -p 
2 (2i/ + 1) (i/ + n) 3 2 I 2z/ + 2, 1 



+ 



z/ + n 



3-^2 



-n, p+1, — p 
2z/, 1 
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which complements relation (3.12) of 

p (p + 1) F ( 1 - n , -P, P + l 
n + 2u 3 2 y 2v + 1, 2 

_ g (p + 1) / 1 - n, 1-p, p 
2u + 1 3 2 ^ 2z/ + 2, 2 

reproduced here for completeness. One needs to derive transformations A3. 15H — (13.170 directly from 
the advanced theory of generalized hypergeometric functions [3], [7]. 

It is worth noting, in conclusion, that explicit solutions of the systems of the first order difference 
equations with variable coefficients are not widely available in mathematical literature. This is why, 
it is important to investigate in detail a remarkable structure of the expectation values pointed out 
in this paper for a classical problem of the relativistic quantum mechanics, such as spectra of high-Z 
hydrogenlike ions. 

Acknowledgments. I thank Carlos Castillo-Chavez, Hal Smith and Vladimir Zakharov for valu- 
able discussions and encouragement. 

Appendix A. Laguerre and Dual Hahn Polynomials 
The Laguerre polynomials are j7], [16], [TT] 



T(a + m + l] 



-m 



Kn (*) = \ v( ' 1F1 ; x . (A.l) 

ml 1 (a + 1) yct + i 



The dual Hahn polynomials are given by [16 J 

«S?(.<.+ 1),«,1) = (1 + a - fc) - ( , 1 + a + C) " (A.2) 




X 3-^2 



1 + a — b, 1 + a + c 



In fl3.ip - fl3.3p and (I3.4p - (j3.6p of this paper, we are dealing only with the following special cases: 
m = n, n — 1 and a = b = 0, c = 2u, s = p and a — b — 0, c — 2i/H, s— p + l, respectively. 

The difference equation for the dual Hahn polynomials has the form 

where Af (s) =Vf(s + l) = f(s + l)-f(s),x(s) = s(s + l), x x (s) = x (s + 1/2) , and 

a(s) = (s - a) (s + b) (s - c) , (A.4) 
a (s) + r (s) Vxi (s) = o" (-s - 1) (A.5) 
= (a + s + 1) (6 - s - 1) (c + s + 1) , 
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A m = m. (A. 6) 

It can be rewritten as the three-term recurrence relation 

a (-s - 1) Vx (s) y (s + 1) + a (s) Ax (s) y (s - 1) (A.7) 
+ (A m Ax 0) Vx (s) Vxi (s) - a (s - 1) Vx (s) - cx (s) Ax (s)) y (s) = 0. 
See jTT], [H] and [16] for more details on the properties of the dual Hahn polynomials. 

Appendix B. Matrix Identity 

The required matrix identity 

/ a (2k + e (p + 1)) + 2/i a (2ek + p + 1) + 2e/x \ 

/i + an pi + a/t / 

a(2re + £(p + 1)) - 2/i a (2ek + p + 1) - 2e/j, ^ 

\ /i — an /i — an / 

-p (4v 2 e + 2k (p + 1) + (2ke +p + l)) 4/i 2 (p + 1) + p (2ks + p) (2ne + p + 1) 
-p (4z/ 2 + 2/te: (2p + 1) + e 2 p (p + 1)) 4/i 2 e (p + 1) + p (2k£ + p) (2k + e (p + 1)) 

(/i + an) (a (2en + p) — 2efi) (/i — an) (a (2en + p) + 2e/z) 
(/i + an) (a (2k, + ep) — 2/i) (/i — an) (a (2k, + ep) + 2/i) 

-a 2 p 2 + 2ae/j,p — 2 (/j 2 — a 2 K 2 ) 2 (/i 2 — cl 2 k 2 ) 

,2 ^2,„2 



B.l) 



x 



X 



= 8a V (p + 1) 

provided a 2 = 1 — e 2 and fi 2 = k 2 — u 2 , can be verified with the help of a computer algebra system. 



-2 (/a 2 — a 2 /t 2 ) a 2 p 2 + 2ae\xp + 2 (/i 2 - 
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